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ABSTRACT 

t;;-!- ' The strong magnetic fields {B ~ 10^^ — 10^'^ G) characteristic of neutron stars make 

O . all the properties of an atom strongly dependent on the transverse component K± 

of its generalized momentum. In particular, the photoionization process is modified 



O I substantially: (i) threshold energies are decreased as compared with those for an atom 

^ ■ at rest, (ii) cross section values are changed significantly, and (iii) selection rules valid for 

atoms at rest are violated by the motion so that new photoionization channels become 



allowed. To calculate the photoionization cross sections, we, for the first time, employ 
exact numerical treatment of both initial and final atomic states. This enables us to take 



;_i ■ into account the quasi-bound (autoionizing) atomic states as well as coupling of different 

c/3 . ionization channels. We extend the previous consideration, restricted to the so-called 

' centered states corresponding to relatively small values of K±, to arbitrary states of 

atomic motion. We fold the cross sections with the thermal distribution of atoms over 
^ ■ K. For typical temperatures of neutron star atmospheres, the averaged cross sections 

differ substantially from those of atoms at rest. In particular, the photoionization edges 
are strongly broadened by the thermal motion of atoms; this "magnetic broadening" 
exceeds the usual Doppler broadening by orders of magnitude. The decentered states 
of the atoms give rise to the low-energy component of the photoionization cross section. 
This new component grows significantly with increasing temperature above 10^'^ K and 
decreasing density below 1 g cm"'^, i. e., for the conditions expected in atmospheres of 
middle-aged neutron stars. 

Subject headings: atomic processes: photoionization, autoionization — magnetic fields 
— stars: neutron 
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1. INTRODUCTION 

The thermal-like surface radiation detected re- 
cently from a number of isolated neutron stars (NSs) 
(see, e. g., Ogelman 1995 for a review) provides 
valuable information about the NS temperatures and 
radii. This enables one to trace the thermal history 
of NSs and to put constraints on the properties of 
the superdense matter in NS interiors (e. g., Pethick 
1992). However, the observed spectra are substan- 
tially modified by reprocessing of the thermal radia- 
tion in NS atmospheres (Romani 1987). Pavlov et al. 
(1995; and references therein) have shown that the 
strong magnetic fields of NSs, B ^ 10^^ — 10^^ G, af- 
fect the emergent spectra significantly. In particular, 
photoabsorption by the strongly magnetized atomic 
hydrogen gas which presumably covers the NS sur- 
face should be taken into account. The magnetic field 
strongly affects the structure and radiative transitions 
of the hydrogen atom when the electron cyclotron 
energy hujc ~ fieB / {nicc) is much greater than the 
Rydberg energy Ry = mce^/{2h'^) = 13.6 eV, i. e., 
7 = hL0c/{2Ry) = 5/(2.35 x 10^ G) > 1. In partic- 
ular, the atom becomes stretched along the magnetic 
field, the ionization potential is increased by a factor 
of ~ In 7, and radiative transitions become strongly 
dependent on the photon polarization. 

Photoionization of strongly magnetized hydrogen 
atoms at rest has been considered by a number of au- 
thors. The most thorough consideration in the adia- 
batic approximation, applicable for very large 7 and 
oj 'ti iOc {^ is the photon frequency), was presented 
by Potekhin & Pavlov (1993; hereafter PP93). This 
approximation neglects coupling of states with differ- 
ent values of the Landau quantum number n which 
specifies the energy of the electron motion transverse 
to the magnetic field (the other quantum number as- 
sociated with the transverse motion, the projection s 
of the angular momentum onto the field, is exact for 
atoms at rest). Potekhin, Pavlov, & Ventura (1996; 
hereafter PPV97) went beyond the adiabatic approx- 
imation and resolved, in particular, a longstanding 
contradiction: some authors (e. g., Schmitt et al. 
1981; Wunner et al. 1983; Ventura et al. 1992) con- 
cluded that the photoionization cross section is ex- 
actly zero for photons polarized transversely to the 
magnetic field, whereas other authors (e. g., Gnedin, 
Pavlov, & Tsygan 1973; Miller & Neuhauser 1991; 
PP93) obtained finite cross sections. PP93 showed 
that the discrepancy rooted in different forms of the 



interaction potential (so-called "velocity form" and 
"length form") employed by the two groups of au- 
thors, and PPV97 proved that the velocity form used 
by the former group resulted in missing the main con- 
tribution to the cross section. The nonadiabatic ap- 
proach enabled PPV97 to study correlation of dif- 
ferent ionization channels and manifestation of au- 
toionization (Beutler-Fano) resonances in the pho- 
toabsorption spectra. 

The above-cited papers did not take into account 
the atomic motion inevitable in any real situation. 
If there were no strong magnetic field, the motion 
would only lead to the trivial Doppler broadening 
of the spectral lines and photoionization edges. On 
the contrary, the finiteness of the atomic velocities in 
strong magnetic fields leads to many qualitatively new 
effects (Pavlov & Meszaros 1993; hereafter PM93) as- 
sociated with strong distortion of the atomic structure 
due the Lorentz electric field generated by the motion. 
Since the moving atom loses its cylindrical symme- 
try, the projection s of the angular momentum is no 
longer a good quantum number (in other words, the s- 
adiabatic approach which neglects coupling of states 
with different s becomes invalid). Hence, the selec- 
tion rules for the radiative transitions arc changed so 
that new types of transitions become allowed. Since 
the atomic structure depends on the velocity of the 
atom, and different atomic states are affected by the 
motion differently, a new type of broadening of spec- 
tral features, magnetic broadening, becomes the most 
important. 

The energies and wave functions of the moving 
atom depend on the transverse component K± of a 
generalized momentum K (see below), which, con- 
trary to the canonical and kinetic momenta, is con- 
served for a neutral atom moving in a magnetic field. 
PM93 analyzed various effects of the motion of the 
hydrogen atom in the perturbative regime, when K± 
is small. When K± exceeds some critical value, the 
magnetic and (motion-generated) electric fields form 
a potential well at a distance Tc = mcK±/LOc from 
the proton, which gives rise to so-called decentered 
(or shifted) states (Burkova et al. 1976). The decen- 
tered states recover the cylindrical symmetry (s be- 
comes a good quantum number again at large K±), 
but the axis of symmetry, with respect to which s is 
defined, is shifted by re- In other words, these states 
can be described in the shifted s-adiabatic approxi- 
mation. Binding energies, wavefunctions and oscil- 
lator strengths of radiative transitions between the 



discrete levels of an arbitrarily moving H atom (be- 
yond the adiabatic approximation) have been studied 
by Potekhin (1994; hereafter P94). Our recent study 
of the hydrogen bound-bound opacity in NS atmo- 
spheres (Pavlov & Potekhin 1995; hereafter PP95) 
was based on these results. 

Less investigated are effects of the motion across 
the field on the bound-free radiative transitions of 
atoms. PM93 have shown, in the perturbation ap- 
proach, that the motion strongly affects photoion- 
ization cross sections through opening and enhance- 
ment of additional photoionization channels forbid- 
den for the non-moving atom due to its cylindrical 
symmetry, and through the strong magnetic broad- 
ening of the photoionization edges. Bezchastnov & 
Potekhin (1994; hereafter BP94) suggested a modifi- 
cation of the adiabatic approach which allows for dif- 
ferent displacements of the initial and final states and 
presented a convenient formalism for calculating the 
transition rates between such states. They have also 
performed calculations with the initial state treated 
non-adiabatically; the two approaches exhibited a fair 
agreement with each other. Kopidakis, Ventura & 
Herold (1996; hereafter KVH96) also went beyond 
the perturbation approximation and calculated the 
frequency dependences of the cross sections for vari- 
ous values of Kj_ , assuming these values small enough 
to neglect the decentered bound states. The results of 
the two latter papers differ from each other by orders 
of magnitude for radiation polarized transversely to 
the field. The origin of the discrepancy is the same 
as for atoms at rest — KVH96 employed the velocity 
form of the interaction potential which led to wrong 
results when used with the n-adiabatic approxima- 
tion which treats the Landau number n as the exact 
quantum number. Both BP94 and KVH96 treated 
final (continuum) states of electron in the (shifted) 
adiabatic approximation, neglecting coupling of fi- 
nal states with different s and contribution of quasi- 
bound (autoionizing) states to the photoionization 
cross section. 

The general non-adiabatic numerical approach de- 
veloped by PPV97 for the continuum wave functions 
of the strongly magnetized hydrogen has been applied 
there to the particular case of non-moving atoms. In 
the present work, we employ that approach, together 
with the theoretical formalism of P94 and BP94, for 
an accurate treatment of both discrete and continuum 
states of the moving atom and the radiative transi- 
tions between them. Unlike BP94 and KVH96, we do 



not restrict ourselves to the case of relatively small 
generalized momenta of atoms, but consider also the 
decentered states which arise at large K±. This en- 
ables us to calculate the cross sections of photoion- 
ization by arbitrarily polarized photons of the hy- 
drogen atom arbitrarily moving in a magnetic field 
B ^ 10^^ — 10^'^ G. Averaging these cross sections 
with thermal distributions of atoms yields the bound- 
free opacities needed for modeling NS atmospheres. 

In Section 2 we outline basic techniques for numer- 
ical treatment of a strongly magnetized H atom and 
its interactions with radiation (the general form of the 
interaction matrix elements is given in the Appendix) . 
In Section 3 we present the photoionization cross sec- 
tions of H atoms in various states of motion across the 
magnetic field and cross sections averaged over ther- 
mal motion for conditions typical for atmospheres of 
isolated NSs. 

2. THEORETICAL FRAMEWORK 

2.1. Generalized momentum and decentering 
of H atoms 

The center-of-mass motion of an atom in a mag- 
netic field can be conveniently described by the gener- 
alized momentum K (Gorkov & Dzyaloshinskii 1968; 
see also Johnson, Hirschfelder, & Yang 1983), 
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where P is the canonical momentum. To — Vp is the 
relative coordinate, the subscript "e" or "p" indi- 
cates electron or proton, respectively, and the cylin- 
dric gauge of the vector potential, A{r) = ^B x r, is 
used. We shall consider the representation in which 
all three components of K have definite values. Then 
the two-body wavefunction can be factorized into a 
phase factor depending on the center-of-mass coordi- 
nate and the wavefunction ipK depending on ro — rp. 
It is useful to define a basic electron-proton sep- 
aration rjVc and to regard the deviation from it, 
r = To — Tp — rjrc, as an independent variable. Here 
Tc = {c/eB'^)B X K IS the position of the relative 
guiding center, and rj is the shift parameter. It is also 
convenient to use a special axial gauge of the vector 
potential (Vincke, Le Dourneuf, & Baye 1992; P94). 
Then we arrive at a set of Hamiltonians ff*-'" and 
wavefunctions ip^'^'{r) depending on ry. In the partic- 
ular cases of ?7 = and rj — I, the "conventional" and 
"shifted" representations of Gorkov & Dzyaloshinskii 



(1968) are recovered. Wavefunctions with different rj 
are related to eacli ottier by the phase transformation 
(BP94) 
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where rj_ = (x, y) is the coordinate transverse to the 
field B, and toh = rnc + rrip is the total mass of the 
atom. 

If there were no Coulomb attraction, then the 
electron Landau number n = 0,1,2,... and the z- 
projection of the angular momentum of the relative 
motion s > —n would be exact quantum numbers. In 
this case the energy of the transverse motion (with 
the constant minimum energy subtracted) is 



E^^ ^ n+ (mo/mp)(n - 
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(3) 



and the transverse part of the wave function can be 
described (for a given K ±) by a Landau function 



<i>„,(rl) = -^ e--^ /„+,,„ U^ (4) 



which depends on r^ — re± — rp± — rc {= r±{rj = 1)). 
Here (j) is the azimuthal angle of r'^, am = {fic/eBY^'^ 
is the magnetic length, /„„/ is a Laguerre function 
(Sokolov & Tcrnov 1968), and the z-projection s of 
the angular momentum of the relative motion is de- 
fined in the shifted reference frame, 77 = 1 (since K± 
is definite, the electron and proton do not possess def- 
inite z-projections of the angular momenta separately 
from each other — see Johnson et al. 1983). Equation 
(0) implies that after photoionization of a moving hy- 
drogen atom, when the photoclcctron moves away to 
large distances at which the Coulomb interaction can 
be neglected, it acquires the displacement Tc from the 
proton in the transverse plane. 

2.2. Atomic wavefunctions 

A wavefunction -0^ of an atomic state \k) can be 
expanded over the complete set of the Landau func- 
tions 

V'i'')(r) = 5]<i>„,(r^)5i'?l;.(z) (5) 

ns 

(note that r± depends on ?/). The adiabatic approxi- 
mation would correspond to retaining only one term 
in this expansion. 



A bound state of the atom can be numbered as 
\i) = |ni. Si, V, K), where rii and Si relate to the lead- 
ing term of the expansion (ph , and v enumerates lon- 
gitudinal energy levels 



E. - E., - E„ 



(6) 



and controls the 2-parity: gn,s;K.{—z) = (— l)''g„'^i;„(z). 
For a moving atom, this way of numbering is unam- 
biguous at 7^300 (P94). 

Since the transverse factors $„s in equation (^ are 

known analytically, only the one-dimensional longitu- 

(77) 
dinal functions gnAn are to be found numerically. For 

the bound states, numerical calculations (Vincke et 
al. 1992; P94) yield the following. At smaU trans- 
verse generalized momenta K±, the states v = Q 
remain tightly bound and centered, with the elec- 
tron cloud concentrated around the proton. For such 
states, the conventional representation (?/ = 0) is ap- 
propriate. At large K^, the states are decentered, 
with the electron localized in a potential well shifted 
apart from the proton. Then 77 = 1 is the apt choice. 
With growing K±, the transition from the centered 
tightly-bound states to the dcccntered states occurs 
in a narrow range of K± near K^ — ^/2m-H_\E^, and 
the shifted (77 = 1) adiabatic approximation becomes 
fairly good at K± ^ Kc- For the hydrogen- like states 
7^ > 1, however, the mean electron-proton separa- 
tion grows steadily, being close to re at arbitrarily 
small K± , so that these states can be described by the 
shifted adiabatic approximattion at any K± . Finally, 
at very large K± (3> ^Ti/a^, where ae is the Bohr 
radius) the longitudinal functions become oscillator- 
like for any states, corresponding to a wide, shallow 
parabolic potential well of a depth ~ e^/rc below the 
continuum boundary (Burkova et al. 1976; P94). 

A continuum state can be numbered as |/) = 
\nf, Sf,Ef, /, K'), where Uf and s/ correspond to an 
open channel for the energy Ef, and / = ±1 corre- 
sponds to a type of the final state — e. g., to a given 
z-parity or to a direction of electron escape, depend- 
ing on asymptotic conditions. In the photoionization 
process the value of K is constrained by the momen- 
tum conservation law: K' = K + hq, where q is a 
photon wavevector. 

For the continuum without motion {K'j_ = 0), a 
technique for calculating the longitudinal wavefunc- 
tions has been developed by PPV97. It can be easily 
extended to the general case, K\ 7^ 0, by numerat- 
ing the orbitals ns in equation (@) for the state |/) 



with single integers j such that larger values of j cor- 
respond to higher energies E^^. Since the Coulomb 
interaction ceases to distort the transverse wavefunc- 
tions at infinity, the orbitals in equation (0) become 
uncoupled at large z provided that the full-shift rep- 
resentation [t]' = 1) of the continuum state is chosen. 
We shall adopt this choice hereafter. 

2.3. Interaction vifith radiation 

The cross section for ionizing an atomic state \i) 
into a continuum state |/) by a photon with frequency 
w, wavevector q, and polarization vector e is 
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where /z is the reduced mass, Lz is the z-extension of 
the periodicity volume of the final state, and M is the 
dimensionless interaction operator. For an arbitrarily 
moving atom of finite mass, the expression for M is 
given in the Appendix. Using equation (0) for both 
the initial and final atomic states, we obtain 



(/|M|z)= J2 {n'^WJl 
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(8) 
where \nsr], k)h corresponds to the longitudinal wave- 
function (7„s;«;(z), and \nsri)±^ corresponds to $„s- 
The interaction operator M in this equation depends 
on 77', ry, because the common phase factor (eq. M) 
has been included in it (see the Appendix). 

The Landau functions for the initial and final states 
in equation (||) depend on differently displaced argu- 
ments if 77' 7^ 77. This is the case for the tightly bound 
initial states at K± < Kc, for which the apt choice of 
the shift parameter is 77 = (while 77' = 1). The dif- 
ferent displacements strongly complicate explicit in- 
tegration in the inner (transverse) matrix elements. 

In order to reduce the initial and final wavefunc- 
tions in the same basis, KVH96 expressed the dis- 
placed Landau state as a series over the undisplaced 
states. This reexpansion involves an additional sum- 
mation index which runs over a huge number of val- 
ues if the relative displacement is large enough, which 
renders the reexpansion method impractical at high 

The difficulty can be circumvented using the creation- 
annihilation operator formalism of BP94. Analyti- 
cal calculation of {n's'ri'\M\nsri)± is performed in the 



Appendix. As a result, the interaction matrix element 
{f\M\i) reduces to a sum of one-dimensional quadra- 
tures (eq. [|A9|) feasible for numerical evaluation. 



3. RESULTS AND DISCUSSION 
3.1. Photoionization at fixed Kj^ 

Many properties of the photoionization cross sec- 
tion can be understood from the Xj^-dependence of 
the level energies shown in Figure H^. The mono- 
tonic increase of the energies shifts the photoioniza- 
tion thresholds redward. Admixture of many orbitals 
near "anticrossings" of levels belonging to different 
s-manifolds (e. g., at K± « 150 a.u.) causes pecu- 
liarities in the radiative transition rates (P94). The 
levels which lie above zero in Figure |l| correspond to 
quasi-bound states in the continuum. They give rise 
to sharp resonances in the cross sections at the pho- 
ton energies corresponding to transitions between the 
initial bound state and the quasi-bound states. 

Figure || demonstrates a comparison of our nu- 
merical results (solid lines) with those obtained us- 
ing approximations encountered in literature. The 
short-dash-dot curves are obtained using the velocity 
form of the interaction operator with the s'-coupling 
"switched-off" (s'-adiabatic approximation). We see 
that this coupling is very important for evaluating 
the cross section (j_ (for the left circular polarization 
transverse to the magnetic field) and less important 
for the other two basic cross sections, cr_|_ (right circu- 
lar polarization transverse to the field) and crii (linear 
polarization along the field). The dotted lines arc ob- 
tained using the shifted adiabatic approximation for 
the final state and the length representation of the 
interaction operator. Note that they are close to the 
short-dash-dot lines in the upper panels. This is in ac- 
cordance with approximate equivalence of using the 
7i-adiabatic approximation (tt., n' = 0) with the length 
form and using the velocity form with allowance for 
n,n' > 0. This equivalence has been suggested by 
PP93 and proved by PPV97 under the condition that 
the non-adiabatic corrections were small. It fails for 
K± ^ Kc if the s'-coupling is neglected (compare 
dotted and dot-short-dash curves in the lower pan- 
els), but we have checked numerically that it holds if 
the s, s'-coupling is included. 

The dot-long-dash curves correspond to the ap- 
proach of BP94. It differs from the plain shifted 
(r7's')-adiabatic approximation (dotted curves) by em- 
ploying the orthogonalization of the final state with 



respect to the initial one (see eq. [32] of BP94). At rel- 
atively small K± (upper panels), this approximation 
yields qualitatively correct results. With increasing 
K±, the approximate (T+ becomes progressively less 
accurate due to violation of the cylindrical symmetry 
of the continuum states of the atom; for instance, it 
exceeds our cross section by a factor of ~ 5 — 7 at 
K± = 100 a.u. However, the approximate cr_ and ctii 
remain surprisingly accurate. 

Short dashes correspond to the approach of KVH96, 
who use the n-adiabatic approximation for the initial 
state, and the shifted adiabatic and Born approxi- 
mations for the final state. If to abandon the Born 
approximation (retaining the adiabatic one), then a± 
would become considerably lower, by 0.5 to 3 orders 
of magnitude for the parameters chosen in Fig. y, 
and even more different from the true cross sections. 
This lowering is due to the approximate orthogonal- 
ity of the initial and final longitudinal wave functions, 
which is not provided by the Born approximation. 

All the discussed approximations miss the Beutler- 
Fano resonances (narrow peaks and dips near the 
thresholds), which occur at the energies of the quasi- 
bound states in the continuum. These resonances 
are analogous to those discussed by PPV97, except 
that now these quasi-bound states all belong to the 
{n' = 0)-manifold. Autoionization in this case be- 
comes possible due to the motion-induced coupling of 
different s'-channels. 

Except for the missing resonances, the discussed 
approximations appear to be sufficiently accurate for 
the longitudinal (crii ) polarization. 

Figure ^ demonstrates a representative sample of 
our numerical results for 7 = 1000. In accordance 
with Figure 1, the photoionization thresholds shift 
redward with increasing Kj_. Since Kc ~ 150 a.u. 
at this field strength, the left and right panels corre- 
spond to the centered and decentered atoms, respec- 
tively. 

For the centered atoms, the slopes of the curves 
for each given polarization remain practically coinci- 
dent with those at K± = 0. The values of the cross 
sections for the parallel and right polarizations, cor- 
responding to the radiative transitions with s' = s 
and s' = s + 1 allowed at K± = (when s is a 
"good" quantum number), decrease only slightly with 
increasing K± from to /"fc- For the left polarization, 
radiative transitions from the ground state are forbid- 
den {a- = 0) at /sTj. = (PP93). With growing Kj_, 



the coupling between different s-channels arises, and 
cr_ grows (approximately as Kj_ — see PM93) until 
it becomes comparable to cr+ at K± ~ Kc- 

For the decentered atoms (right panels of Fig. W, 
the energy dependences steepen with growing K±. 
This is explained by the change of the effective poten- 
tial from Coulomb-like at K± < K^ to oscillator-like 
as the atom becomes decentered. The effective poten- 
tials corresponding to different s become nearly iden- 
tical, that results in a significant decrease of cr+. Still 
faster is the decrease of a- caused by restoring ap- 
plicability of the s-adiabatic approximation and cor- 
responding selection rules at large K±. At the same 
time, the threshold value of (T|| remains almost inde- 
pendent of Kj_ . 

Since the moving atom loses its cylindrical symme- 
try, the cross sections depend on the angle if between 
K± and the transverse component of the wavevec- 
tor. The (^-dependence is most pronounced when 
both the wavevector q and polarization direction e 
are perpendicular to B. The cross sections cr|| and a± 
are drawn in Figure ^ by dot-long-dash and dashed 
lines, respectively. The former cross section does not 
depend on ip in the dipole approximation, whereas 
the small non-dipole corrections (eq. AlOc |) cause a 
weak dependence with an amplitude <C 1%. There- 
fore, we show only one curve for an on each panel. 
For the transverse polarization, the (^-dependence be- 
comes pronounced at K± ^ Ji/am ~ 30 a.u. The up- 
per and lower dashed curves correspond to the largest 
and smallest possible values of aj_{ip), at (p — and 
(fi = 90°, respectively. The reason is that the trans- 
versely polarized photons "see" the largest transverse 
size of the atom (stretched perpendicular to K±) 
when they propagate along K±, and the smallest 
one when propagating transversely to K±. The (p- 
dependence between these two extremes, illustrated 
in Figure 0, is described by the formula 
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Note that in the absence of coupling of different s- 
channels the amplitude A in equation (H) would equal 
^cr+(T_, which is always greater than the actual value 
of A, calculated numerically. 

The curve corresponding to Ej — 1.1 Ry in Fig- 
ure goes higher than that corresponding to Ej — \ 
Ry because the ionization channel s' = \ opens at 
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hui 
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1.089 Ry, where uj, 



{■mc/mp)uc 
is the proton cyclotron frequency. The curve corre- 



cp 



spending to Ef — 2 Ry goes much higher than the 
other curves because it faUs onto a slope of a Beutler- 
Fano resonance. Such resonances are most prominent 
at K± ~ Kc (cf . Fig. 0) because of the violation of the 
s-adiabatic approximation in this range of the trans- 
verse generalized momenta. Nevertheless, these peaks 
remain very narrow. Figure ra presents some autoion- 
ization widths Fa calculated according to equations 
(14) and (15) of PPV97 (with obvious generalization 
to the case of moving atoms) . The left panel shows 
the widths of the tightly-bound states s = 1 through 
4 which enter the continuum at K± ^ 1200, 560, 340, 
and 237 a.u., respectively (see Fig. ||). Owing to the 
similarity of the effective potentials related to differ- 
ent s values, the jumps of the curves corresponding 
to opening of additional autoionization channels occur 
at nearly the same values of K± . The hydrogen- like 
levels of the (s = l)-manifold (right panel of Fig. ^) 
belong to the continuum at any K±. The autoion- 
ization widths vanish at K± ^ K^ and at K± ^ K^ 
because of restoring cylindrical symmetry of the wave- 
functions. 

3.2. Absorption by atoms in thermal equilib- 



The bound-free absorption coefficient for a hydro- 
gen gas of given temperature and density is propor- 
tional to the photoionization cross section CTk^k' {K±) 
folded with Nf^{K), the distribution of the number 
density of atoms in the initial state. Since the level 
energy Ef^{Kj_) tends to a constant value at K^ -^ oo, 
the normalization integral for the Boltzmann distri- 
bution N^{K) ex exp[-EjKj_)/kBT] diverges. The 
divergence is eliminated if nonideality of the gas is 
taken into account. At high K±, the binding energies 
become small, and atomic sizes become large. This 
means that the atom can be easily destroyed by sur- 
rounding particles, so that very high K± should not 
contribute to the corresponding integrals. In order 
to allow for this effect, we use the occupation prob- 
ability formalism, described in detail by Hummer & 
Mihalas (1988) and generalized to the case of the K±- 
dependent atomic structure by PP95. Following this 
approach, we write 



N^{K) (X wM<i_)exp[-E^{K^)/kBTl 



(10) 



where Wk{K±) is the occupation probability, equal to 
the fraction of atoms (in the specified state \k,K)) 
which essentially preserve quantum-mechanical (par- 
ticularly, optical) properties of an isolated atom. We 



estimate 'w^{K±) according to equation (14) of PP95, 
based on the quantum-mechanical atomic sizes cal- 
culated by P94. A rapid decrease of w^{K±) at 
K± > 'jfi/aB (see Fig. 3 of PP95) limits contribu- 
tion from atoms with large generalized momenta and 
provides convergence of integrals over K±. 

Since the threshold energies decrease with increas- 
ing K±, folding the photoionization cross sections 
with Nii{K±) leads to the redward broadening of the 
photoionization edges, with a typical width T^ ~ 
ksT (PM93). It has been shown by PM93 and PP95 
that the "magnetic width" is much greater than the 
Doppler and (impact) coUisional widths under typical 
conditions of NS atmospheres. In this paper we take 
full account of the magnetic broadening and neglect 
other broadening mechanisms. 

The averaged photoionization cross sections for 
atoms on the ground level are shown in Fig. |6| for two 
values of temperature T. For comparison, we show 
also cross sections for the bound-bound (b-b) tran- 
sitions from the ground level (PP95). Both the b-b 
and bound-free (b-f) averaged cross sections are cor- 
rected for the stimulated emission, i. e., multiplied by 
(1 — exp[— Tioj/fesT]). For the circular polarizations, 
the b-b and b-f transitions dominate in different (ad- 
jacent) spectral ranges — the low-energy b-f compo- 
nent, arising from the decentered states (cf. Fig. ^, 
turns out to be much weaker than the b-b component 
in the same energy range. For ctm , this is also the 
case at logT = 5.2, but not at logT = 5.8. In the 
latter case, the low-energy b-f component dominates 
the absorption spectrum, except at very low energies 
hoj ^ 40 eV. 

The peak of a+ in the right panel of Figure o at 
huj « 15 eV and the corresponding peaks of the 'b-b' 
curves arise near the proton cyclotron frequency Wcp, 
around which the transitions of strongly decentered 
atoms are concentrated (cf. PP95). 

Figure R demonstrates the effect of density on 
the averaged photoionization cross sections. With 
increasing p, the low-energy component rapidly de- 
creases because the decentered atoms, due to their 
large sizes, are easier destroyed by the plasma mi- 
crofields. The growth of the high-energy component 
noticeable at higher T is explained by the relative 
increase of the number of centered atoms when the 
decentered atoms are destroyed. It should be noted 
that at high densities one more effect, broadening of 
the photoionization edges by quasi-static interactions 
with surrounding particles, may become important. 



As the density grows, highly excited bound states 
become randomly disturbed by the surrounding and 
form a quasi-continuum. This results in continuous 
photon absorption below the undisturbed photoion- 
ization threshold (e. g., Dappen, Anderson, & Miha- 
las 1987). These "bound-quasi- free" transitions re- 
quire a separate treatment and are not included in 
the present paper. 

Figure shows the effects of T and B on the 
photoabsorption spectral shapes. With growing T, 
the low-energy components significantly increase, and 
the high-energy components of a\\ and cr+ decrease 
slightly, because the decentered states become more 
populated. At lower T, the high-energy component 
of CT_ grows with T linearly. 



TujJcv, ksT 



-'cp 

^^000,010 



(7+ 



(11) 



in accordance with the result of PM93 obtained in 
the perturbation regime (eooo,oio is the energy of the 
main b-b transition for the right polarization - e. g., 
£000,010 ^ 70 eV for 7 = 1000). At higher T, the 
growth of a^ decelerates because of decreasing frac- 
tion of the centered states and nonperturbative ef- 
fects. 

At higher B, the above-discussed proton cyclotron 
resonance for a+ becomes stronger. Owing to the 
approximate similarity of the energy spectra related 
to different s-manifolds (see Fig. |^), the individual 
Beutler-Fano resonances, although smoothed by the 
thermal averaging, are collected together into regu- 
larly placed peaks. Such peaks corresponding to four 
proton cyclotron harmonics are indicated by arrows 
for (T- at 7 = 3000 in Fig. R At high densities, the 
peaks can be further smoothed by interaction with 
surrounding particles. 



mal modes for any frequency and direction of propa- 
gation should be found for solving the radiative trans- 
fer equations. Bulik & Pavlov (1996) have described 
how these polarizations can be calculated if the basic 
cross sections a± and a\\ are known. 

Our results show that atomic motion in strong 
magnetic fields drastically changes both the bound- 
bound and bound-free opacities. In particular, it con- 
siderably broadens the spectral lines and photoioniza- 
tion edges (with a typical magnetic width ~ ksT), 
gives rise to additional absorption at low frequencies 
due to decentered atoms, and opens new transition 
channels forbidden for atoms at rest. These unusual 
properties of the spectral opacities alter the structure 
of the NS atmospheres and manifest themselves in the 
spectra and light curves of radiation emergent from 
the NS surface layers. 

The authors are grateful to Victor Bezchastnov 
and Joseph Ventura for fruitful discussions, and to 
Cole Miller and the referee, Tomek Bulik, for use- 
ful remarks. This work was partially supported 
through NASA grant NAG5-2807, INTAS grant 94- 
3834, RBRF grant 96-02-16870a, and DFG-RBRF 
grant 96-02-00177G. AYP gratefully acknowledges 
the hospitality of Joseph Ventura at the University 
of Crete, Gilles Chabrier at ENS-Lyon, and Chris 
Pethick at Nordita, where a part of this work has 
been done. 



4. CONCLUSIONS 



The results presented here provide, together with 
the results of PP95 on the bound-bound transitions, a 
basis for calculations of realistic opacities of NS atmo- 
spheres. To find the opacities, one should add tran- 
sitions from excited atomic levels, which is straight- 
forward. More difficult are the problems of ioniza- 
tion equilibrium in the dense, non-ideal plasmas and 
closely related problem of the non-ideality effects on 
the bound-free transitions (including transitions to 
the quasi-free states below the unperturbed photoion- 
ization threshold). Finally, polarizations of the nor- 



APPENDIX: INTERACTION MATRIX ELEMENTS 

Velocity and length representations of the operator of interaction of a moving hydrogen atom with radiation 
have been derived by BP94. Here we employ the velocity representation as the simpler one. KVH96 have noted 
that the interaction of the radiation field with the electron spin should be taken into account. It has been found 
by PPV97 that the spin term almost does not change the cross sections at energies below 10 keV. Nevertheless, we 
include this term for the sake of generality, neglecting the terms responsible for the spin-flip transitions (Schmitt 
et al. 1981) which are strictly forbidden at the considered photon energies below hujc- 

In the conventional representation of the wavcfunctions, equation (A6) of BP94 yields for the dimensionless 
interaction operator M: 



(2aB)"^M = cxp i—^qr 

' TOH 



rric K q 



:(gxe) 



Top II 



771c K 

TO-H h 



e, 



(Al) 



where q is the photon wavevector, e is its polarization vector, and negligibly small terms ^ {mc/mp)q are omitted. 
Here r = Vg — rp is the relative coordinate, and tv and II are the electron- and proton-type kinetic momentum 
operators: 



P 



2c 



-B X r, 



n = n B X r. 

^ 2c 



Their cyclic components act on the Landau states as 



■!T±\ns)± 
n±|ns)_L 



n* |n± l,s =F 1)_L, 



T — 
am 

ih 

T — Vn* + s\n,sT l)j_, 



(A2a) 
(A2b) 

(A3a) 
(A3b) 



where n» = n+l and n^ = n for the right (+) and left (— ) components, respectively. If the full-shift representation, 
rj' = 1, is used for the final state, and an arbitrary shift parameter rj for the initial one, then, according to equation 
(A15)ofBP94, 



M 



2h 

— exp 



1 771, 

2~ 



"7H 



Ki-7?)js:rj 



e (exp [i(qr_L/2 + g^z)] F^ + exp \-vf]rcq - iqr^/2] Fp) 



(A4) 



where r ~ r^^ r^ — r\rc- Once again, the terms ~ qmc/mp are omitted. The components of the vector operators 



Fe and Fp are 



Fc+ = — + (1 - V) + , 

?77c 7MH fnc 



Fp±- 
F, 



n± 

TOp 
1 



(1 - 7;) 



K± 



-{Pz 



TOH 

%./2), 



7r_ , ,K_ 
i^c- = — + (1 - r?) , 

TOc TOH 



Fpz ^ Pz/mp. 



(A5a) 
(A5b) 

(A5c) 



In order to proceed with equation (|^), we need to calculate the transverse matrix elements {n's'r]'\M\nsri)±. In 
the particular case of the dipole approximation {q — > 0), such calculation has been performed by BP94. In the 
general case, we use the following result obtained with the same operator technique as in BP94: 



{n's'r]\Kj 



q±|exp 



1 7(7, 

2~ 



TOH 



— {v' - v)Kr±_ + iqr^/2 



\nS'q,K\_)±_ = Jnn' iCc)Jn+s,n'+s' i^p), (A6) 



where q is arbitrary, 



Co = -iflm ( ^^iv' - v)K+/h + iq+ + rj'q+)/2 



ma 



and 



(A7a) 
(A7b) 

(A8) 



Jn„'(|?|e''^)=e'("-"')^/„„,(ien, 
Inn' being a Laguerre function (Sokolov & Ternov 1968). 

Using equations (A4), (A3a), (|A3bD, and ([Aq) for rj' — 1, and neglecting an insignificant common phase, we 



finally arrive at the expression for the interaction matrix element (^: 



{fim = 2^j2 K^^n' 



s'ns ^ ^-^^^n's'ns ^ ^ z ^V-l- „i g' ns 



(A9) 



where 



AI^V] 



[\'nJn-l,n'{^c,+ ) + Cc,-Jnn'{£,c, + )\ Jn+s,n' +s' [^p,+ ) ^n' s' 






exp(-i77rcq) [Vn + S + 1 Jn+s+l.n'+s'{^p-) ~ Cp-Jn+sm' + s'i^p-)] 



(p) 



X Jnn'i^c,-) ^n's'ns^ 



^'^n's'ns — [Vn+1 Jn+l^n'{£,c, + ) - Cc,+ Jnn'{(,c, + )] Jn+s,7i' +s' {^p,+ ) Z,^, ^, 



(c) 



m, 



+ —2- exp(-i77rcq) [\/n + S Jn+s+l,n'+s' (?p,-) + Cp. + Jn+s,n'+s'i^p 



X '^nn'yS.c.-J ^n's'ns^ 






r(<=) 



r(c) 



^n's'ns ' n^rn'dz^n's'r 



+ exp(-i77rcq) 



TOn 



^n7i' vSc,+ /<-'n+s,n'+s' \.sp, + j 



?(P) _ la o Z^P^ 

^n's'ns 2 ^^^ "'*'"*i 



In these equations 



Cc,+ 

Cp,± 



7i \mH 



7i TOfj 



{1^V)K-, 



h toh 



(l-^)^±, 



^e,± and ^p^± are given by equations (A7) with rj' — 1 and q — ±q, respectively, and the factors 



?:('=) - /«'»' f\Ji-^ 



(p) 






(AlOa) 

(AlOb) 
(AlOc) 

(Alia) 
(Allb) 

(A12a) 
(A12b) 



are the longitudinal matrix elements subject to numerical evaluation. 

We have presented the result for rj' = I. For arbitrary 77', the matrix elements are obtained in analytic form 
using equation (A15) of BP94 and our equations (A6), (AT) with q — [2mc + ri'{mp — mc)]q/mH and q = — [2mp + 
r]'{me — 7TT-p)]q/mH. The result is a straightforward generalization of equations (AlO). 
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